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STUDY OF A CONTROL SYSTEM TO ALLEVIATE AIRCRAFT 
RESPONSE TO HORIZONTAL AND VERTICAL GUSTS 

By William H. Phillips 
Langley Research Center 

SUMMARY 

An analysis is made of the longitudinal response of an airplane in three degrees of 
freedom, with provision for including the effects of a vane -controlled gust-alleviation 
system and a simple pitch stability augmentation system. The systems are designed to 
compensate for the effects of vertical and horizontal gusts. Examples are presented to 
show the effect of such systems with various combinations of parameters on the response 
of a STOL airplane with externally blown flaps in the landing-approach condition. The re- 
sults include transient responses after penetration of step horizontal and vertical gusts, 
responses to step elevator deflections, roots of the characteristic equations, and magni- 
tudes of the gust forcing terms. 

The complete elimination of gust forcing terms in the longitudinal-response equations 
of an airplane, which is possible when horizontal gusts and speed changes are neglected, 
does not appear possible with systems of the type studied when these factors are taken 
into account. Nevertheless, the important forcing terms and the resulting response may 
be reduced to relatively small values by use of flaps with suitable aerodynamic character- 
istics, a rearward center- of-gravity location, and gearing of flaps to spoilers or other 
drag devices to reduce longitudinal force changes due to flap deflection. The stability 
characteristics of the alleviated airplane appear to be unsatisfactory without stability 
augmentation to provide pitch -attitude stabilization and improved damping of the phugoid 
mode. 


INTRODUCTION 

A previous analytical study of methods of reducing the vertical acceleration due to 
gusts (ref. 1) was limited to consideration of vertical gust inputs, and the airplane dynam- 
ics were simplified by assuming constant airspeed. These restrictions are justifiable for 
high-speed flight conditions. For STOL airplanes, particularly in the landing- approach 
condition, the relative sensitivity to horizontal gusts is increased, and effects of speed 
variations may be more critical because of the reduced period of the phugoid motion. For 
this reason, an extension of the theory to include speed variations and horizontal gusts is 
desirable. 



The gust-alleviation system of reference 1 utilized a combination of wing flaps and 
elevators geared to move together so that, for analysis purposes, they could be considered 
as one control. These surfaces were moved in response to an angle- of- attack vane mounted 
ahead of the nose in such a way as to offset the forces and moments applied to the airplane 
by the gusts. In this sense, the system operates as an open- loop control. The vane position, 
however, is also affected by airplane motion. The resulting feedback of airplane motion to 
the control surfaces affects the longitudinal stability. The objective of the study of refer- 
ence 1, therefore, was to reduce the gust inputs while retaining adequate stability and con- 
trol characteristics. This objective was accomplished in the cruise flight regime without 
employing additional stability -augmentation devices. One objective of the present study is 
to investigate the possibility of offsetting the effects of both horizontal and vertical gusts, 
taking into account speed variations as well as vertical and pitch responses. The response 
of a STOL airplane in the landing- approach condition is used as an example to investigate 
this system. 

The method of analysis used in reference 1 was based on the conventional linearized 
airplane stability theory, in which a first-order approximation to lag effects is employed. 
This method is advantageous because the addition of a gust- alleviation system does not in- 
crease the order of the equations above that of the basic airplane equations, and because 
most of the effects of the system may be visualized as modifications to the stability deriva- 
tives of the basic airplane. The modified derivatives provide physical insight into the ef- 
fects of the alleviation system. Relations were derived in reference 1 which must be satis- 
fied by the design parameters in order to eliminate completely the response to vertical 
gusts. The usefulness of these relations was later demonstrated in a research flight pro- 
gram on a gust- alleviated airplane (ref. 2). Extension of the theory to include horizontal 
gusts is desirable to determine whether similar relations can be found which will permit 
reduction of the response to horizontal gusts. 

In the present investigation, therefore, longitudinal equations of motion for the air- 
plane are set up in three degrees of freedom, considering both horizontal and vertical 
gust disturbances. A simple gust-alleviation system is assumed, consisting of a spring- 
loaded vane-type sensor operating wing flaps through a servomechanism. The spring- 
loaded vane-type sensor when used with a gust -alleviation system has the capability of 
providing gust alleviation over a range of speeds without the need for gain changes between 
the sensor motion and flap response. The vane -type sensor, although it is useful for 
research purposes, has practical disadvantages compared with internal instruments, such 
as accelerometers and gyros, because it is relatively fragile and exposed to damage. The 
relations derived for this type of system, however, are useful in providing a guide for the 
design characteristics of systems using internal sensors. 

The analysis presented herein neglects flexibility effects and is based on the as- 
sumption of a conventional configuration in which forces and moments on the stabilizing 
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surfaces can be calculated separately from those on the wing-fuselage combination. More 
detailed analysis would be required to investigate effects of structural flexibility and con- 
trol-system dynamics in an actual configuration. 

SYMBOLS 


A aspect ratio 


normal acceleration, g units 
L 


lift coefficient, 


v 2 

■&- S 


C L equilibrium lift coefficient 


C pitching-moment coefficient, 
m 7 


M 


pV 


Sc 


'mow win f?~kody pitching-moment coefficient at zero lift with flaps deflected in 
’ ’ approach condition 


Cx longitudinal-force coefficient, 


X 


PV 


Ci vertical-force coefficient, 




C z equilibrium vertical-force coefficient 
c mean aerodynamic chord of wing 


d 

D differential operator, — 

ds 

F ratio of elevator deflection to pitch angle 
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ratio of elevator deflection to Dtf 


ratio of flap deflection to vane deflection 
ratio of radius of gyration to chord, k y / c 

radius of gyration in pitch 
lift 

ratio of tail length to chord 

ratio of distance between vane and center of gravity to chord 
pitching moment 
spring moment on vane 
mass of airplane 

ratiq of flap deflection to elevator deflection with vane fixed 
period, sec 

nondimensional pitching velocity, 0c/2V 
wing area 

V 

distance traveled in chord lengths, — t 

c 

time for oscillation to damp to one-half amplitude, sec 

time for oscillation to double amplitude, sec 
time 


increment of horizontal velocity 



u 


11 


u g 

v 

w 

*; 

X 

z 


nondimensional horizontal velocity increment, u’/V 
horizontal gust velocity 

nondimensional horizontal gust velocity, u^/V 
true velocity with respect to undisturbed air mass 
weight 

vertical gust velocity 

increment in vertical velocity with respect to undisturbed air mass 
longitudinal force, positive forward 
vertical force, positive downward 
increment in angle of attack 

1 

w 

increment in angle of attack with respect to undisturbed air mass, — - or d - y 

V 


a equilibrium angle of attack, - y 

y increment in flight-path angle 

y equilibrium flight -path angle 

S e increment in elevator angle from equilibrium position 

Sj increment in flap deflection from equilibrium position 

8 v increment in vane angle from equilibrium position 

S variation of vane angle with horizontal velocity, 

v u ^ 

d s 

8 variation of vane angle with angle of attack, — I 

v a 3a 
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e increment in downwash angle 

l damping ratio 

6 increment in pitch angle 

6 equilibrium pitch angle 

f± relative density factor, m/pSc 

p air density 

r' time constant of servomechanism 

r nondimensional time constant of servomechanism, X r', chords; or time 
constant of aperiodic mode, chords c 

nondimensional time constant of servomechanism for pitch stability augmen- 
tation system, chords 

w frequency, rad/ chord 

Subscripts : 

b value for basic airplane 

g gust 

o with respect to inertial axes fixed in undisturbed air mass 

t tail 

v vane 

w wing-fuselage combination 

Dot over quantity indicates differentiation with respect to time. 
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Stability derivatives are indicated by subscript notation; for example, 



dZ_ 

d a 

W 


C 


z 


a 


W 


SC Z 

3a 

W 


3C r 


C Z 


Du 


uc 

2V 



'^Z 

0C 

B 2V 


Subscript after a stability derivative indicates component of airplane which contrib- 
utes the derivative; for example, ^C Zg 


THEORETICAL ANALYSIS 

The relative importance of vertical and horizontal gusts in producing changes in 
wing lift may be obtained from the relation 


AL = 


“ t t 

u 

-iC r + -i (2 C t ) 


V 


'a V 




Thus, the ratio of the incremental lift for horizontal and vertical gusts is 


AL H 2c l 


aL v c l 


The mean- square values of the horizontal and vertical gust velocities are expected to be 
about equal. The ratio ALjj/ALy in cruising flight may therefore be in the range from 
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0.1 to 0.25. In the landing-approach condition, however, the value may exceed 1.0. Design 
of gust- alleviation systems for STOL airplanes should therefore include consideration of 
horizontal gusts. 


Description of System 

A sketch of the system under consideration is given in figure 1. The spring-loaded 
vane, which is discussed in detail subsequently, is sensitive to both horizontal and vertical 
gusts. The output of this vane operates the flaps through a servomechanism. The flaps 
are assumed to be capable of upward and downward motions from their initial deflected 
position. If necessary, gearing other controls to the flaps to modify their characteristics 
is considered within the scope of the analysis. Such controls might include the elevators, 
spoilers, drag brakes, or engine nozzles. In fact, one of these types of controls, such as 
spoilers, used alone might prove satisfactory. In this case, the same analysis can be 
applied by substituting the symbols and aerodynamic characteristics appropriate to spoiler 
control in place of those for flap control. The pilot’s control input is also fed into the flap 
control system to provide direct lift control, inasmuch as the alleviation system prevents 
generation of lift for maneuvering through changes in angle of attack. 


Equations of Motion 


The longitudinal equations of motion are written first in dimensional form, consider- 
ing separately the contributions of the wing and tail. The equations of motion are derived 
with respect to wind axes. The wing flaps are assumed to be available as controls. The 
axis system and angles used are given in figure 2. The equations are as follows: 


mu' - u X , + u!X , + a + O..X + S »X j. +&X„ 

O W L W t cij. I 5 | o 


m <*0 - V»> - "w Z U^ + “t Z uf + “w Z . w + * t Z a t + 8 f Z S f * « Z f 


mky6> = u!..M, 


uIM.. 


w-u^ + ^-uj. 


a,„M +a,.M + 8 .M t 

W a t a ^ I e 


+ 9 M, 

f 


To nondimensionalize the equations, the following substitutions are made: 
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D _ d _ c d 
ds V dt 


fi - 


m 

PSc 


C X = 


pV 


c z 


£Z!s 


M 


m 


V 2 

Sc 


The equations then become 


2mDu q -u w C + u t C X + a w C X + a t C X + s f C X, +0C 'i 

Ui a a, b 4 


W 


w 


f 


w 


M D ~ “ D9 j ~ u w C Z + U t C Z + a w C Z + a t C Z + S f C Z c . +eC Z- 
\ V / u w u t a w a t o f e 


2 uK D 0 = u C + u, C + a C + ci, C + 8 ,*C + 8 C 

Y W m U t m u, W m a * m a + f m 8, m 0 

w t w t f 

In the subsequent analysis, all variables are expressed in terms of u Q , u g , a Q , 
and 8 and their derivatives. Lag or lead in the occurrence of a velocity increment 
with respect to the motion at the wing is expressed by use of Laplace transform notation. 
The required relations are as follows: 


u = u„ + u„ 

W 0 g 


w' 


w’ 


tt w" a o + a g 


u, = u + u„e 
to 


-ID 


cie -fD - ID --&D ~ D p p.,, 

o. , — a - a e + a„e -a — e - 8 . e x + l DO 

t 0 0 3a S g3a f 
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A vane sensor which is sensitive to the vertical and horizontal increments of velocity 
at the vane location is assumed to be mounted ahead of the nose. The velocities at the vane 
are given by the expressions 

ID 

u = u + u e n 

V o g 

ID 

a v =a o +a g e ~ D0l n 

The vane is assumed to respond to these velocity increments with no lag. A control law 
is now introduced which relates the motion of the wing flaps to the motion of the vane, with 
a constant time lag approximating the lag of the servomechanism which actuates the flaps: 


8 f = KS e" rD 
f v 

where 

S v = 5 v U v +S v “v 
u a 

It will be shown subsequently that for the type of vane considered, s = -1. Hence 

v a 

S f = K(s u -a )e" rD 
f v u V V 7 

The lag or lead terms i, and r are expressed in chords. These quantities may be 
converted to dimensional quantities in seconds by multiplying them by the ratio c/V. 

In the subsequent examples, the elevators, spoilers, or drag devices may be as- 
sumed to be geared directly to the wing flaps. For purposes of analysis, the effect of 
these devices may be taken as a modification of the flap derivatives, inasmuch as their 
motion is in phase with that of the flaps. 

The values for u y and may be substituted in the expression for Sp which is 
in turn substituted in the expression for a fc . All expressions for lag or lead are then 
represented by the linearized approximation: 

rD i 

e = 1 + rD 
or 
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Discussion of the effect of this approximation is given in reference 1. The expressions 
for u w , u^, a w> a. ^ ; and § f are substituted in the equations of motion. After combining 
terms multiplying each variable or its derivatives, the equations of motion may be written 
in the following form: 


2,aDu - uC v - i DuC y - a. C„ - I D a C v -0 C Y - - D8C„ 
u 2 A Du ° Xa o 2 0 X Dc 0 X 6 2 X q 


~ D 2 £C Y 9 = a C +1 D a C + u C Y + i Du C v 

4 X D 2 0 ^ X 2 g X Da g X u 2 g 


(la) 


g 


g 


g 


2^ (Da - D0) - uC 7 - - DuCr, - a C„ - - Da C 7 ~ec„ - - D£ 
° 2 ^Du o o 


D0C r 


a o 2 u "Da o 2 “q 


1 r.2. 


- TT9C 7 p - a C + ± Da C ry + U C„ + _ Du c 7 
4 Z D 2 6> B z a _ 2 8 Z Da B Z u 2 8 Z Du 


(lb) 


g 


g 


g 


2mK. 2 D 2 6> - uC - 1 DuC„ - a_C_ 


y 


m. 


u 2 Du 


o m 


-i- Da C -6C - - D&C 


“o 2 


° m Da m e 2 


rri 


Id 2 ^C 9 = a_C_ + — D a C + u C + I Du C (lc) 

4 m D*d B m a _ 2 B m Da g m u 2 g m Du 


g 


g 


g 


The C z derivatives are given as follows: 


Be 


C 7 = Cy =C +Cy -KS ^1C 7 + KS Cy 

" \ “w \ V “ 3S f “t V “ Z *f 


(2a) 


pZn =KS v If U + T)C Z ' KS v " C Z S 

2 DU U U 


(2b) 


B € ,xr Be 


C Z =C Z = C z + C„ / 1 - + K 

l a 

O g 


KC r 


a - a - a w a t \ 3a BS. 


(2c) 


* c 

2 Z Da 


0 L 


l — - Kii (r+ l) 

Ba 3S 


C Z + KtC z 


J a t 


(2d) 
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1 

2 



U - Ki 


n 



K ' t n C Z. 


(2e) 


1 

4 


C z 2 = JL 
Z D 4? n as 


(r + -t)C r 


- Kri C„ 
n Zr c 


(2f) 


1 

2 





— <*n 

38 n 


r - i) C z - K(^ n - t)C z 

J a t *f 

(2g) 


1 

2 







•«Z * K 5 v<V'> C Z t 
a t U S f 


(2h) 


The C v and C m derivatives can be obtained by substituting X and M for Z 

ill 

wherever it occurs. 


The 9 derivatives come from the change in the components of gravity forces acting 
on the airplane resolved along the wind axes resulting from an increment in 9, The 
expressions for these derivatives are as follows: 


X Q = W COS y 

C x = -C L cos y 

6 


Zq = W sin y 



- -C^ sin y 


Effect of Feedback of Pitch Angle and Pitch Rate to the Elevator 

A simple stability augmentation system sensitive to pitch angle and pitch rate is 
considered in the subsequent analysis. The effect of this system may be analyzed, to the 
same degree of approximation as that used in the preceding derivation, by assuming the 
control law: 


S e = (F 0 + GD>9) (1 - r 0 D) 
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The nondimensional gain G relating S 0 to D6> may be converted to a dimensional gain 
G’ relating & e to 8 by the formula 


G’ = G _ 

V 

Likewise, the lag t q expressed in chords may be converted to a lag in seconds by 
the formula 


r e T e 


The increments to the C z stability derivatives provided by this stability augmentation 
system are as follows: 


’ F % 


el t 


AC„ = 2 (G-tF)/c 


J D o 


8 'f Zg 


e/t 


aC z d 2 0 = ' 4Gt * c z 8 


e/t 


Again, the increments to the C x and C m derivatives can be obtained by substituting 
X and m for Z wherever it occurs in these expressions. 


Vane Characteristics and Control- System Parameters 
Required for Gust Alleviation 

The forcing terms for the complete three degrees of freedom which are associated 
with gust effects may be collected from equations (1) as follows: 


a. C-y + D <1 ( C-y - Cy ) + U Cy + Du Cy 

8 x a 2 S X Da X g Xu 2 8 X Du 

& 


(3a) 


a C17 + — Da fc 17 - C»7 \ + U Cy + — Dll C7 

B z a 2 8 Z Do Z q s z u 2 8 Z Du 0 


(3b) 
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(3c) 


a C 
g m. 


+ 1 Da 


m 


- C 


Da 


m. 


+ u C + I Du C 
« m u 2 6 m Du 


g 


The response to vertical and horizontal gusts may be eliminated by reducing the 
coefficients of all these terms to zero. It is desirable, therefore, to determine how many 
of these terms may be reduced to zero by use of the system under consideration. 

The response to vertical gusts under conditions of constant airspeed may be elimi- 
nated if the single-underlined terms are reduced to zero. This condition may be provided 
if the following relations, derived in reference 1, are satisfied: 


K = 




m < 


m 


w 


K 



3 K 


(4a) 


(4b) 


(4c) 


This result may be verified by substituting these relations in the a and Da terms of 

& § 

expressions (3) and using the values of the derivatives given in equations (2). 

The required value of C m may be obtained with little performance penalty by 

8 f 

gearing the elevator to move in phase with the flaps. The required value of _L, however, 

3S f 

is opposite in sign from that normally associated with wing flaps. In the research flight 


program described in reference 2, the reversed sign of 



was obtained by gearing an 


inboard segment of the flaps to move in the opposite direction from the flaps covering the 
remainder of the wing span. This method reduces the overall lift effectiveness of the 
flaps and might be undesirable for a STOL airplane in which flap deflection is required 
for improved maximum lift coefficient as well as for gust alleviation. Nevertheless, the 


a e 



ability to adjust 
14 


to the desired value will be assumed in the present analysis. 



The coefficient of the double-underlined term in expression (3b), C~ , may be 

u 

reduced to zero by use of a vane with the correct value of S in conjunction with the 

v u 

values of K and _1L given in equations (4a) and (4c). By setting the expression for this 

3S f 

coefficient (eq. (2a)) equal to zero, the required value of s is found to be 

v u 


( Cz u)b 

( Cz »)b (5) 

where the values of (c r ^ j and ^ re P resen t values for the basic airplane before 

the addition of the gust-alleviation system. 

The required value of 8 may be obtained by use of a suitable spring moment on 

v u 

the vane, tending to make it float down. The vane will then move up in response to either a 
rearward horizontal or an upward vertical gust. Consider the vane subject to velocity and 
angle- of -attack variations. The spring moment on the vane is equated to the aerodynamic 
moments by the relation 


M c = (8 + a)C 

s v m 


v 2 

tL. s c 

S 2 v v 

V 


whence 


= 


m < 


U 


pv* 


S c 

V V 


a 


The spring moment is assumed not to vary with vane deflection. 
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The variations of vane angle with angle of attack and with velocity are then 



Equation (6) may be solved for M as follows: 

S 


( 6 ) 


M 


s 



S c 

V V 


This relation establishes the required strength of the spring moment on the vane, tending 
to make it float down. 


The preceding formulas were derived on the assumption that the variation of lift with 
velocity is approximated by a linear relation for small changes in speed from the trim 
condition. This type of analysis might be required for a STOL airplane, for which the 


power effects on 


( Cz u). 


are large. If the lift at a given angle of attack varies as v 


as it does in the case of a rigid airplane with negligible slipstream effects, then the value 


of 


^ e Q ua l s 2 C z , where is the equilibrium value of the vez'tical-force 


coefficient. The value of the spring moment is then 


C Z C m 


= 


S v 2 

Mb 


S c 
v v 


- 2 

Inasmuch as C 2 varies inversely as V , the spring moment required for correct al- 
leviation of vertical and horizontal gusts is then constant throughout the speed range. 

If one additional condition is specified, the triple-underlined terms of expressions 
(3b) and (3c) are also reduced to zero. This condition is 
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The value of C 7 depends on the tail load for trim, which the designer may control by 

N 

varying the value of pitching moment of the wing-flap combination about the moment ref- 
erence point or by varying the center-of-gravity location. Even if all these conditions are 
satisfied, however, the important terms u„C m and all the forcing terms in the longi- 
tudinal force equations remain in the equations and are changed but not eliminated by the 
alleviation system. 


NUMERICAL EXAMPLES 


In order to determine the importance of the forcing terms remaining in the equations 
of motion, a number of examples have been studied utilizing the characteristics of a STOL 
airplane with externally blown flaps in the landing- approach condition. The parameters 
were estimated from reference 3. 


Because the effects of some variables, such as flap deflection, were based on wind- 
tunnel measurements at only two deflections, these parameters are only approximate. The 
values are believed, nevertheless, to incorporate the main characteristics which show 
marked differences from those of a conventional airplane. These differences include, for 


example, large effects of power on the derivatives / C 


and 


u b 


m. 


u/b 


and large 


changes in drag due to flap deflection. The physical characteristics and derivatives for 
the basic airplane are given in table I; The stability and gust-response characteristics of 
the basic airplane and of the airplane with a gust-alleviation system of the type analyzed 
have been studied by calculating the roo'ts of the characteristic equations, the coefficients 
of the forcing terms caused by gust inputs, and the transient responses to step vertical 
and horizontal gust inputs and to pilot control inputs. All calculations were made with a 
high-speed digital computer. These results are presented for six cases, as follows: 
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Case 1: 
Case 2: 
Case 3: 
Case 4: 
Case 5: 


Case 6: 


} 


Basic airplane; c.g. at 0.33c 
Alleviated airplane ; c.g. at 0.33c 
Basic airplane; c.g. at 0.594c 
Alleviated airplane; c.g. at 0.594c 
Alleviated airplane; c.g. at 0.594c; 


} 




= 0 


Alleviated airplane with pitch autopilot; 
c.g. at 0.594c; C x =0 


Response to 
gust inputs 


Figure 3 
Figure 4 


Figure 5 


Response to 
pilot control 
inputs 


Figure 6 


Figure 7 


Figure 8 


For the alleviated airplanes of cases 2, 4, and 5, the derivatives are the effective values 
as determined from equations (2), taking into account the contributions of the alleviation 

system. Values of the gain K, the vane floating tendency 8 , and the lag r are 

v u 


determined from equations (4a) and (5). The values of C m 


and _ d ,.i are selected to 

B S f 


satisfy equations (4b) and (4c). As a result, the single- and double-underlined terms in 
expressions (3b) and (3c) are eliminated, except for the effects of slight numerical errors. 
No attempt has been made to meet the condition of equation (7) which reduces the triple- 
underlined terms of expressions (3b) and (3c) to zero. An exception is made in case 6, in 


which the values of C? and 0 1 are selected to give a small stable value of C 

m s f B 8 f ra , 


of -0.45 rather than zero. The values of these parameters are listed in table II. Values 
of the roots of the characteristic equation for the six cases are given in table HI. Values 
are given for the characteristic time in chords of real roots and for the damping ratio and 
frequency in radians per chord of oscillatory modes. In addition, the values are presented 
as time to one-half or double amplitude and period in seconds for the assumed flight speed 
of 32.61 meters per second (107 feet per second). 


Values of the derivatives in the forcing terms on the right-hand side of the equations 
of motion (eqs. (1)) are given in table IV. 


Transient responses for step inputs of vertical and horizontal gust velocity are given 
in figure 3 for the airplane with and without the gust- alleviation system for a center-of- 
gravity location at 0.33c and in figure 4 for the center-of-gravity location at 0.594c, 

Similar responses for the center-of-gravity location at 0.594c with two types of control 


systems incorporating a value of C, 


of 0 are shown in figure 5. The magnitude of 
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the gust velocity for both vertical and horizontal gusts is 0.570 meter per second (1,87 
feet per second). At the flight velocity of 32.61 meters per second, a vertical gust of this 
magnitude produces a change in flow angle of 1°. The changes in angle of attack a Q and 
velocity V are with respect to inertial axes fixed in the air mass and do not include the 
increment of gust velocity. 

Transient response to step elevator deflections of 1° are shown in figures 6, 7, and 
8 for the same conditions as those for the gust inputs. In the cases with gust-alleviation 
systems, the flaps are operated as a function of the longitudinal control input to provide 
direct lift control in the manner shown in figure 1. The elevator is assumed to respond 
without lag to the pilot's control input, whereas the flap is driven through a servomechanism 
with a lag r used in the gust-response calculations. With the vane fixed, the steady- 
state relation between the flap and elevator deflection is given by 

S, = m' s 

f e 

The value of m’ is taken to be -2.18 in all the six cases studied. The value of b> e in this 
formula is the portion of the elevator deflection resulting from direct linkage to the con- 
trol stick and does not include the portion resulting from linkage to the flap to adjust the 
flap pitching moments. 


RESULTS AND DISCUSSION 

An overall picture of the effects of the systems investigated can be obtained from 
the transient response calculations (figs. 3, 4, and 5). These responses were calculated 
for relatively small gust inputs; however, because of the linearity of the analysis, they 
may be multiplied by the appropriate factor to correspond to a larger input. The results 
for the basic airplane show that with the step gust inputs assumed, the changes in airplane 
acceleration are quite small and the need for an alleviation system for this STOL-type 
airplane in the landing- approach condition might be questioned. In landings made in high 
winds, however, particularly in the wake of buildings or other obstructions, severe turbu- 
lence may be encountered so that the resulting airplane accelerations may be objection- 
able. For this condition, the possibility of exceeding the range of validity of the linearized 
analysis must be kept in mind. The linearized analysis, however, is useful as a first step 
in showing the required characteristics of a gust -alleviation system. 

For the forward center- of- gravity location (fig. 3), the alleviation system decreases 
the initial normal accelerations when the airplane encounters gusts but greatly increases 
the response of the poorly damped phugoid mode. The alleviation system also causes an 
initial response in the negative direction for horizontal gusts. For the rearward center- 
of-gravity location (fig. 4), the alleviation system again decreases the critical normal 
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accelerations. The subsequent motion is slowly divergent. When the alleviation system is 
operating, both vertical and horizontal gusts cause an increase in forward inertial speed. 
This horizontal acceleration is a result of the decrease in drag when the flap is deflected 
upward. The initial horizontal acceleration due to the vertical gust is about 60 percent of 
the initial vertical acceleration of the unalleviated airplane. Such large horizontal accel- 
erations due to gusts would probably be quite objectionable to the passengers. 


In order to reduce the horizontal accelerations, a case was tried, again for the 
rearward center-of-gravity location, with the value of C v reduced to zero (case 5). 

8 f 

This condition could possibly be obtained by gearing the flaps to drag devices or spoilers, 
so that the drag devices extend when the flaps move up. Alternatively, spoilers alone, 
trimmed up in the equilibrium approach condition, might provide the desired changes in 
lift with small changes in drag. The other flap derivatives C 7 and C’ were 

** s , m s r 


assumed to be unaffected when C x 


was reduced to zero. Although in practice these 
f 


derivatives would be affected, an alleviation system equivalent to that assumed could always 
be attained, at least for some range of flap deflections, by modifying suitably the gearing 
between the flaps and elevators and by recalculating the flap gearing and downwash 
parameters. Because aerodynamic data for a suitable drag device were not available, 
the details of such an arrangement were not investigated. As shown in figure 5, the 
responses to the vertical and horizontal gusts for this case are almost eliminated. 


An explanation of the behavior shown by the transient responses may be found by 

examining the roots of the stability equations and the forcing terms. Consider first the 

roots (table m). For the basic airplane with a forward center of gravity (case 1), for 

which the static margin is 0.277c, the normal short period and phugoid oscillation modes 

are encountered. The phugoid oscillation has poor damping and a period somewhat 

shorter than expected for this flight speed. The positive value of C tends to shorten 

m u 

the period. For the basic airplane with a rearward center of gravity (case 3), for which 
the static margin is 0.042c, the short period mode becomes a pair of rapid subsidences 
and the phugoid period increases considerably. 

When the alleviation system is operating, the airplane has one real root with a very 
long time constant; this indicates a condition of near-neutral stability. This neutral 
stability is due to the choice of parameters to eliminate the gust forcing terms 

a g 

and C , which are also equal to the values of C™ and C associated with 
m 9 ^ L ni 

a g a a 
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airplane motion. In addition, when the center of gravity of the airplane is at the rearward 
location, the phugoid mode is divergent and has a very long period. In practice, the pilot 
would probably interpret this motion as a straight divergence in airspeed. 

The amount of disturbance of the airplane, even if it is neutrally stable or unstable, 
depends on the magnitude of the gust forcing terms given in table IV. The values of C z , 

a. 

g 

C , Crr and C are reduced to negligible values by the alleviation system. 

m z n m r> 

a Da 

g g g 

At the forward center of gravity (case 2), the value of C is increased by the gust- 

m 

g 

alleviation system and is responsible for the large disturbance of the phugoid mode en- 
countered in this case. At the rearward center of gravity (case 4), however, the value of 
C m is greatly reduced. The value of this derivative, as shown by equation (2a), depends 

U g 

on the value of C for the basic airplane and on the value of C . The positive 
m u m s f 

value of C for the basic airplane with the forward center- of- gravity location may be 
m u 

explained as follows. The negative pitching moment of the wing with deflected flaps in- 
creases less rapidly than the square of the speed because a large part of the moment 
comes from the slipstream, whereas the balancing tail pitching moment varies as the 
square of the speed. As the center of gravity is moved rearward, the balancing tail load 
required is decreased; this results in a smaller positive value of C m . The value of 

C is also greatly reduced as the center of gravity moves back to a location near the 
“*f 

center of lift of the flaps. In the configuration under consideration, both of these effects 
are reduced to zero at a center of gravity near 0.594c, the rearward center- of- gravity 
location assumed in this study. 


Despite the reduction in the value of C m , other forcing terms remain to produce 

a response to gusts of the alleviated airplane with the rearward center -of- gravity location. 
In view of the increase in forward speed noted previously in this case, it was thought that 
the large positive value of C x was primarily responsible for the remaining disturbance. 

U g 

This value may be reduced by reducing the value of Cy . For the case shown in figure 5, 


with the value of 



arbitrarily reduced to zero, the responses to vertical and hori- 


zontal gusts are almost eliminated. Even in this case, some forcing terms remain in the 
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equations. In particular, the value of C x remains close to that of the basic airplane. 

a 

This derivative, however, apparently has only a small effect on the response. 

Although the gust disturbances on a STOL airplane may be greatly reduced by the 
techniques studied, the stability characteristics without additional stability augmentation 
appear unsatisfactory. The very large time for subsidence of one of the real roots indi- 
cates that the airplane is indifferent as to pitch attitude. In other words, with C and 

near zero, the airplane is in equilibrium at any angle of attack and the lift is inde- 

a 

pendent of angle of attack. Any slight application of elevator control would cause a con- 
stantly changing pitch attitude. 

The response to pilot control inputs is of importance in the overall evaluation of a 
gust -alleviation system. In the alleviation systems studied, a rapid response in normal 
acceleration to a longitudinal control input is provided by the direct-lift control feature 
in which the flap system moves in response to a pilot control input. The subsequent mo- 
tion, however, is affected by the stability characteristics of the alleviated airplane. 

As shown in figure 6, the response of the basic airplane with the forward center-of- 
gravity location (static margin, 0.277c) to an up-elevator control deflection is weak, with 
the normal acceleration reversing after a few seconds because of the decrease in airspeed. 
A somewhat greater initial response is obtained with the alleviation system and direct- 
lift control operating, but the poor damping of the phugoid mode may lead to control 
difficulties. 

The response of the basic airplane with the rearward center-of- gravity location 
(static margin, 0.042c) appears unduly sensitive because of the low static stability (fig. 7). 
The airplane continues to diverge upward in pitch after the normal acceleration reverses. 
Normally, an airplane of this type would not be flown with such rearward center of gravity 
without stability augmentation. 

The rearward center-of-gravity location was tried, in this case, to obtain improved 
gust alleviation. As shown previously, the airplane with a large drag increment due to 
flap deflection experienced an undesirable divergence in airspeed due to a gust disturbance. 
The longitudinal control characteristics appear undesirable because of the unusually large 
instability of the phugoid mode. Because of the instability and the long period, this motion 
would probably be interpreted by the pilot as a straight divergence. 

The alleviated airplane with C x =0 (fig. 8) has somewhat less objectionable 

s f 

characteristics. The phugoid mode is again unstable, but the rate of divergence is reduced. 
Because of the effectively zero values of C m and C for the alleviated airplane, the 

m ct “u 
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airplane continues to pitch up at a fairly constant rate as long as the elevator control is 
applied. Another case, for which data are not presented, was studied with a value of C m 

a 

of -0.45 which corresponds to the value of C of the basic airplane at the rearward 

center-of-gravity location. The results were not appreciably changed. 

The foregoing results indicate that although the gust disturbances on a STOL airplane 
may be greatly reduced by the techniques studied, the stability characteristics without 
additional stability augmentation appear unsatisfactory. In order to investigate the effects 
of stability augmentation, a simple system sensitive to pitch angle and pitch rate was 
added to the alleviated airplane with C-^ - 0. The sum of these feedback quantities was 

S f 

added to the motion of the elevator, as indicated in figure 1. The effects of this type of 
control on the gust response are shown in figure 5 and on the response to a control input 
in figure 8. 

The response to gust inputs is still essentially eliminated by the alleviation system 
with the pitch stability augmentation system (fig. 5). The response to a control input ap- 
pears to be more desirable for this case than for the other cases studied (fig. 8). The 
pitch angle responds in the correct direction but does not increase at a constant rate and 
the reduction of airspeed is more gradual. The roots (table III) are all aperiodic conver- 
gences except for the near-zero real root, which is a very slow divergence. 

In all the examples presented, increased normal acceleration is accompanied by a 
fairly rapid decrease in airspeed. This characteristic is typical of airplanes operating at 
high values of lift coefficient. The pilot would be required to coordinate throttle action 
correctly with elevator control to make a correction to the flight path. The application 
of control decoupling techniques to the gust-alleviation system studied would appear to 
provide a means of further improving the longitudinal control characteristics. 

As stated previously, an internal sensor, such as an accelerometer, has practical 
advantages over a vane. The flap motions on encountering a gust with such a sensor are 
similar to those with the vane. It therefore appears desirable to employ flaps having 
aerodynamic characteristics which reduce the gust forcing function just as in the case of 
the vane sensor. The use of flaps having these characteristics was shown in reference 1 
to be desirable for a transport-type airplane in the cruise condition with an accelerom- 
eter sensor. Such flaps, in addition to producing a change in lift, should have pitching- 
moment and downwash characteristics similar to those shown to be required with the vane 
sensor. In addition, the use of a rearward center-of-gravity location to reduce C m and 

C and the gearing of spoilers or drag brakes to the flaps to reduce would be 

ra s f 8 f 

desirable. 
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A rearward center-of-gravity location requiring stability augmentation to provide 
adequate longitudinal stability also appears advantageous from the performance stand- 
point. The down load on the tail required for trim, which is very large for STOL airplanes 
with large trailing-edge flaps, is thereby reduced or reversed; this results in the require- 
ment for less tail area as well as less wing area and thrust for the same landing-approach 
speed. 


CONCLUDING REMARKS 

This study considers the alleviation of response to horizontal and vertical gusts of a 
STOL airplane utilizing a vane-type sensor which operates flaps and other controls through 
an automatic control system. The response to horizontal gusts is shown to be important 
for a STOL airplane in the landing-approach condition. 

The complete elimination of gust forcing terms in the longitudinal response equa- 
tions of an airplane, which is possible when horizontal gusts and speed changes are ne- 
glected, does not appear possible when these factors are taken into account. Nevertheless, 
the important forcing terms may be reduced to relatively small values by use of suitable 
flap aerodynamic characteristics, a rearward center-of-gravity location, and gearing of 
flaps to spoilers or drag devices to reduce longitudinal-force changes due to flap deflection. 

The stability characteristics of the alleviated airplane appear unsatisfactory without 
the use of stability augmentation to provide pitch-attitude stabilization and improved 
damping of the phugoid motion. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 21, 1973. 
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TABLE I.- AIRPLANE CHARACTERISTICS AND DERIVATIVES 


USED IN ANALYSIS 
(a) Physical characteristics 


Dimensional data: 
W 


m 

c 

S 

St 

h 

V 

p 


Nondimensional parameters: 

l 

^n 

A 

M - . 

K y 

C L 

^m,o,w 

9e/95f 

9e/9ai 


245444 N (55178 lbf) 

.... 25028 kg (1715 slugs) 

3.203 m (10.51 ft) 

74.45 m 2 (801.36 ft 2 ) 

33.17 m 2 (357.0 ft 2 ) 

3.667 m (12.03 ft) 

. . . 32.61 m/sec (107 ft/sec) 
1.226 kg/m 3 (0.002378 slug/ft 3 ) 

3.50 

4.09 

7.22 

85.63 

1.144 

5.058 

-2.32 

0.202 

0.430 


Aerodynamic center of wing-fuselage combination 


0.213c 



TABLE I.- AIRPLANE CHARACTERISTICS AND DERIVATIVES 
USED IN ANALYSIS - Continued 

(b) Contributions to stability derivatives of wing-fuselage combination and of tail 
( c Xq.)w 0.920 

( C Z a )w - 8 - 02 

(%)t "O ’ 494 

( C Z„)t -l- 72 

( C X u )w -°- 120 

( C Zu)w - 5 - 84 

( c X u ) t -0.620 

( C Z u )t 0.330 

( CX ^w - 1 ' 92 

( CZ ^ W - 4 ' 30 

(about aerodynamic center of wing-fuselage combination) -1.83 

\ °V w 
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TABLE I.- AIRPLANE CHARACTERISTICS AND DERIVATIVES 
USED IN ANALYSIS - Concluded 


C X« • • 

C Za * * 

Cm <x‘ • 
l r 

2 X D a 

| Cz D a 

ir 

2 tm Da 


Cx d 

Cz e 

Cm e 



C X U 

C Z 


u 


c m u • • 
! C XDu 
? C ZDu 
| Cm Du 


(c) Stability derivatives of basic airplane 


c.g. at 0.33c 
0.638 
-9.00 
-2.49 
-0.743 

-2.59 

-9.06 

-5.06 

0 

0 

-1.73 

- 6.02 


- 21.1 

-0.740 

-5.51 

1.86 

0 


0 


0 


c.g. at 0.594c 
0.856 
-9.00 
-0.376 
-0.169 

-2.59 

-9.06 

-5.06 

0 

0 

-0.393 

- 6.02 

- 21.1 

-0.023 

-5.64 

-0.0809 

0 

0 

0 
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TABLE II.- VALUES OF GAINS AND PARAMETERS FOR ALLEVIATED AIRPLANES 


■ 

Case 

c.g. 

K 

8 v 

u 

r 

m’ 

° m s 

S f 

36 
3 8 f 

2 

0.33c 

1.86 

-0.612 

4.09 

-2.18 

-1.34 

-0.306 

4 

0.594c 

1.86 

-0.627 

4.09 

-2.18 

-0.202 

-0.306 

5 

0.594c 

1.86 

-0.627 

4.09 

-2.18 

-0.202 

-0.306 

a 6 

0.594c 

1.89 

-0.612 

4.09 

-2.18 

0.0386 

-0.266 


cl 

In addition, for case 6, the parameters of the pitch stability augmentation system 
are: F = 1.00; G = 3.41; r Q = 0.50. 


TABLE m,- ROOTS OF CHARACTERISTIC EQUATION FOR CASES STUDIED 


Case 

£ and co, rad/chord, for 
oscillatory modes orr, 
chords, for aperiodic modes 

Tj/ 2, sec, and P, sec, for 
oscillatory modes or Tj/ 2 or T 2 , 
sec, for aperiodic modes 

1: Basic airplane; 
c.g. at 0.33c 
(static margin, 0.277c) 

1 1 = 0,774 co^ - 0.119 

1 2 = 0.034 « 2 = 0.040 

T l/2 » 0.733 p = 5.20 

T 1/2 = 48.6 P = 15.5 

2: Alleviated airplane; 
c.g. at 0.33c 

r 1 = 4270 t 2 = 3.41 

£ = -0.046 w = 0.0307 

T l/2 = 296 T l/2 = 0,234 

T 2 - 48.2 p = 20.2 

3: Basic airplane; 
c.g. at 0.594c 
(static margin, 0.042c) 

= 22.9 r 2 = 7.24 

£ = 0.089 = 0.013 

T l/2 = 1.56 T l/2 = 0.493 

T 1 / 2 = 58.2 p = 47.6 

4: Alleviated airplane; 
c.g. at 0.594c 

tj_ = 4515 t 2 = 5.24 

£ = -0.76 co= 0.00776 

T^^ 2 — 308 ^1/2 = ^*558 

T 2 = 11.5 P = 79.8 

5: Alleviated airplane; 
c.g. at 0.594c; C x =0 
6 f 

= 8220 t 2 = 5.03 

£ = -0.200 co = 0.010 

T l/2 = 562 T l/2 = 0,344 

T 2 = 32.3 P = 61.8 

6: Alleviated airplane 
with pitch autopilot; 
c.g. at 0.594c; C v =0 
f 

t x = -1497 r 2 = 231 

^2 = 10.55 = 5.28 

T 2 = 102 Tj/g = 15.78 

T l/2 = 0.718 T t y 2 = 0.360 
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TABLE IV.- VALUES OF DERIVATIVES IN FORCING TERMS 

OF EQUATIONS (1) 

(a) Longitudinal- force terms 


Case 

c 

'"X 

a. 

g 

- C 
2 X n 

Da g 

C Y 

A u 

g 

2 %u g 

1: Basic airplane; 

0.638 

0.985 

-0.74 

2.17 

c.g. at 0.33c 
2: Alleviated airplane; 

4.48 

0.00145 

1.61 

1.57 

c.g. at 0.33c 
3; Basic airplane; 

0.856 

0.224 

-0.023 

-0.339 

c.g. at 0.594c 
4: Alleviated airplane; 

4.48 

0.00033 

2.25 

-0.479 

c.g. at 0.594c 
5: Alleviated airplane; 

0.919 

0.00033 

0.0170 

-0.479 

c.g. at 0.594c; 





O 

X 

cy? 

Hh 

11 

o 






(b) Vertical 

-force terms 



Case 

C Z 

a g 

— c 

2 C Z n 
JDa^ 
g 

n 

c z 

u g 

i 

1 r 

2 z % 

1: Basic airplane; 

-9.00 

3.43 

-5.51 

-1.15 

c.g. at 0.33c 
2: Alleviated airplane; 

-0.023 

0.0050 

-0.014 

-3.25 

c.g. at 0.33c 
3: Basic airplane; 

-9.00 

3.43 

-5.64 

-0.07 

c.g. at 0.594c 
4: Alleviated airplane; 

-0.023 

0.0050 

-0.014 

-2.84 

c.g. at 0.594c 
5: Alleviated airplane; 

-0.023 

0.0050 

-0.014 

-2.84 

c.g. at 0.594c; 





C Y =0 
S £ 
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TABLE IV.- VALUES OF DERIVATIVES IN FORCING TERMS 


OF EQUATIONS (1) - Concluded 
(c) Pitching- moment terms 


Case 

a s 

X Q 

2 C m 

Da g 


r c m n 

1 : Basic airplane ; 
c.g. at 0.33c 

-2.49 

12.0 

1.86 

-11.6 

2: Alleviated airplane; 
c.g. at 0.33c 

-0.00415 

0.017 

3.38 

-19.0 

3: Basic airplane; 
c.g. at 0.594c 

-0.375 

12.0 

-0.081 

- 2.45 

4: Alleviated airplane; 
c.g. at 0.594c 

0.0001 

0.017 

0.154 

- 9.96 

5: Alleviated airplane; 
c.g. at 0.594c; 

c x . 0 

S f 

0.0001 

0.017 

0.154 

- 9.96 
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Figure 1.- Diagrammatic sketch of control system. 
(Linkages are not drawn to scale.) 








Time , t, sec T' me ’ t, sec 


Figure 3.- Response to step vertical and horizontal gust inputs of 0.570 meter per 
second (1.87 feet per second) for a STOL airplane with and without gust-allevi 
ation system and with c.g. at 0.33c. 
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Figure 5.- Response to step vertical and horizontal gust inputs of 0.570 meter per 
second(1.87 feet per second)for a STOL airplane with two types of gust-allevi- 
ation systems, one of which includes a pitch autopilot. C v = 0: c.g. at 0.594c. 

AS 
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Basic airplane 

— — A1 leviated airplane with direct-lift control 





Time, t, sec 

Figure 7.- Response to a step ele- 
vator deflection of l°for a STOL 
airplane with and without a gust- 
alleviation system and with c.g. 
at 0.594c. 
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Figure 8.- Response to a step elevator deflection 
of 1° for a STOL airplane with two types of 
gust- alleviation systems incorporating direct 
lift ^control, one of which includes a pitch 
autopilot. = 0; c.g. at 0.594c. 
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